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Abstract 

In this work we study boundary crossing probabilities of general one-dimensional diffusion 
processes with respect to constant boundaries. We first give an explicit formula of the prob- 
abilities for diffusions with piecewise linear drift and constant diffusion coefficient in terms 
of its Laplace transform. Then we provide an approximation scheme to obtain the boundary 
crossing probabilities for more general diffusions. Finally, an estimate of convergence rate of 
the approximation is also given. 
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1 Introduction 

Boundary crossing probabilities for stochastic processes play an important role in many research 
fields. For example, in business and industry a company's financial status can be represented by 
its total asset and debt at any given time, which can be regarded as two stochastic processes. 
Then the company goes bankrupt when its asset process falls below its debt process for the first 
time. Therefore the company's default risk can be described by the probability that the difference 
of the asset and debt processes reaches a certain threshold. In finance, the arbitrage-free price 
of a barrier option is given by the expectation of the payoff function of strike price and the first 
hitting time of the underlying asset price to the barrier. In neuroscience, a popular integrate- 
and-fire model assumes that a neuron fires a spike when the membrane potential reaches a 
threshold. In engineering, the system reliability can be represented by the probability that a 
damage process crossing a threshold. In epidemiology, if the spread of an infectious disease is 
modeled by a stochastic process, then it is of interest and importance to know the probability that 
this process reaches a certain threshold within a certain time limit. The first passage time also 
arises in many other disciplines such as biology, chemistry, ecology, economics, environmental 
science, genetics, physics as well as statistics. Some references for applications can be found 
in Wang and Potzelberger (2007) |12| . More recent literature includes Touboul and Faugeras 
(2007) [ID], Taillefumier and Magnasco (2010) [9], Molini et al (2011) [6]. 

Given its important role in many applications, however, the computation of boundary cross- 
ing probabilities is a challenging task. Explicit formulas exist for only a few special instances. 
For example, the reflection principle of one dimensional Brownian motion yields an explicit for- 
mula for it to cross a constant boundary. Further applying the reflection principle in depth leads 
to explicit formulas for linear and piecewise linear boundary crossing probabilities ([3| IT1 fTT]). 
In general it is difficult to obtain explicit formulas for a stochastic process crossing even the 
simplest constant boundary. However, in real applications, the underlying processes are usu- 
ally described by more general diffusion processes than Brownian motion, while the later is 
used to represent random perturbation in a system. Therefore, it is important to study the 
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boundary crossing probability for general diffusion processes. An attempt is made by Wang 
and Potzelberger (2007) [12], who proposed a transformation approach to compute the bound- 
ary crossing probability for a class of diffusions which can be expressed as piecewise monotone 
functional of standard Brownian motion. For this kind of diffusions, explicit formulae for the 
boundary crossing probabilities for certain boundaries can be obtained. However, from both 
theoretical and practical points of view this class of processes is very limited. More specifically, 
|12j consider the following diffusion process 

dX t = ii(t,X t )dt + a(t,X t )dW t , X = x, 

where (Wt)t>o is the standard Brownian motion on R, and show that under some regular 
conditions on fi and a Xt can be expressed as a monotone functional of standard Brownian 
motion if and only if the drift coefficient fi and diffusion coefficient a satisfy the following 
second-order differential equation: 



d r 1 da d /Ida /i 
dxladt dx\2dt a 



0. 



For instance, when a is a constant, n(t,x) must be a linear function of x in order for the above 
differential equation to hold. Hence, the class of processes satisfying the above equation is 
very limited. Moreover, it seems difficult, if at all possible, to generalize this transformation 
method to obtain the boundary crossing probabilities for other diffusion processes that are not 
a functional of Brownian motion. This motivates us to find a novel approach to this problem. 

To simplify notation in this paper we will focus on the diffusion processes satisfying the 
following stochastic differential equation (SDE) with additive noise: 

dX t = n(X t )dt + dW t , X = x. (1.1) 



There is no loss of generality in the sense that any diffusion process satisfying 

dX t = /M(X t )dt + a(X t )dW t 



(1.2) 
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with differentiable and non-zero cr(x) can be transformed into one with unit diffusion coefficient 
through the transformaiton 

Jy o ( ) 

for some i/q. Indeed, it is easy to verify by Ito's formula that the process Yj := F(Xt) satisfies 
SDE HUD with drift coefficient /i(x) = (p/a - a'/2) o F' 1 ^). 

The idea is to first obtain the explicit formula of Laplace transform of boundary crossing 
probability for a diffusion process with piecewise linear drift, and then use this formula to 
approximate the corresponding probability for more general processes. In our approach, an 
estimate of the convergence rate of the approximation can also be derived. More precisely, 
consider the diffusion Xt with continuous and piecewise linear drift fJ,(x) in the form: 

oo 

M0)= Yl 1 [*<-i,*<]( a; )( o i aj + 6 0> (1.3) 

£=— oo 

where sup{|oj| + i = 0, ±1, ±2, • • • } < oo and {. . . < x_ m < . . . < x~\ < xq < x\ < . . . < 
x m < . . .} is a partition of R. We show that the first passage time 

r c = wf{t > 0; X t > c} (1.4) 

admits a density f c (t) on [0, oo) and an explicit formula for its Laplace transform can be obtained 
using the method of [2|. Then for a diffusion process Xt satisfying (II. ID with fi(x) admitting 
bounded derivative, we construct a sequence of processes (X^)^ defined by 

dx\ n) = pn(xi n) )dt + dW t , X^ ] = x (1.5) 

with continuous piecewise linear drifts [i n for n G N such that sup xgR \^{x) — fJ> n ( x )\ < C/n for 

(n) 

some constant C > 0. Finally, we show that X\ converges to Xt almost surely for each t > 0, 

(n) 

and further the boundary crossing probability of X f converges to that of X t in the sense of 
P(sup(xi n) - c) < 0) — ► P(sup(X s - c) < 0), as e -> 0. 

s<t s<t 

In addition, we show that the rate of this convergence is of order 1/n. 
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2 The process with piecewise linear drift 

In this section, we first prove that the boundary crossing density exists for a diffusion with 
piecewise linear drift (jl.3p . Then we give the explicit formula of its Laplace transform using the 
method of Darling and Siegert [2]. 

Given a probability space (£1, F, P), let x € R and (Wf)t>o be a Brownian motion starting 
at Wq = x. For the simplicity of notation, subsequently we will write Wf = W t . Let C x ([0, oo)) 
be the continuous path space starting at wo = x, i.e. 

C x ([0,oo)) = {w : [0, oo) — >• R continuous; wo = xj. 

Let Qa; be the law of Brownian motion (W~t)t>o starting at Wo = x in C x ([0,oo)). Further, 
define F t = cr(w s ; < s < t). By standard theory of stochastic differential equations (eg, [S]), 
there exists a unique non-explosive solution (X t ) t > to the SDE (jl.ip for rt(x) being in the form 
(|1.3p . In order to prove that r c has a density w.r.t. the Lebesgue measure, we shall generalize 
the result of Downes and Borovkov [H Theorem 3.1] to the case where is not necessarily 
everywhere differentiable. 

Lemma 2.1 Let p(s,x,y) denote the transition density of the process (X s ) s >q and q(s,x,y) 
the transition density of the Brownian motion. Let P x denote the law of (X s ) s >q starting at 
Xq = x, and the law of (X s )o< s <j starting at Xq = x and conditioning on Xt = z. Similarly, 
let and Qiji denote the corresponding probability measure with respect to the Brownian motion 
(W s ) s >o. Then for t > and every A £ Ft, 

p*f A ) < ^hlA e G(z)-G( X )-3Mt/2 Q zj A)j M 

p[t,x,z) 

where G(y) = f* Q n(z)&z for some fixed y , and M = inf^gR {/x(y) 2 + with (J,'_(y) = 

liminf 2 ^„^ziiM. 

z—y 

Proof. Since X t = Wt + J n(X s )ds and fj, is bounded, by Girsanov's theorem, 

P*(A) = Q x (CtU), A e F t , (2.2) 
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where Q x {( t lA) = J Ct^-A^Qx and 

Ct = exp [ f fi(X s )dW B - \ f fi(X s ) 2 ds 
L Jo 1 Jo 

Conditioning on the value of corresponding processes at time t, equation (|2.2p can be written as 

/oo roo 
P*(A)p(t,x,z)dz = / P z x (A)P x (X t e dz) 
-oo J ~oo 

/oo 
Q z x (CtU)q{t,x,z) dz. 
-oo 

By the arbitrariness of A £ Tti we have 



K(A) = #^4q^(C,1a). 



p(t,x,z) 



Further, let p(y) = l/y/2%e y2 l 2 and /u n (y) = J 



oo nu(z) 
-oo 



(2.3) 



dz for n £ N. Then n n (y) 



converges to /i(y) for every y £ R as n — > oo. Moreover, by the dominated convergence theorem, 
e(J - ^(X s ))diy s ) 2 = E J (jiniX,) - fi(X s )fds — ► 0, as n ^ oo. 

It follows that, up to an extraction of subsequence, / yu n (X s )dWg — )• / /j,(X s )dW s almost 

Jo Jo 

surely. Similarly by the dominated convergence theorem, G n (X t ) converges almost surely to 
G(Xt) as n — > oo, where G n (y) = J^ q fi n (z)dz. Further, by Fatou's lemma, we have 

Hn{y + h) - n n {y) 



Vn(y) = Jim 



h->0 h 

oo 



> 



lim ~ / (u(y + h--)-u,(y--))-^=e z2/2 d; 

rty-Z + h)-rtv-Z) i - z y 2 



liminf — - — ^—— e - z ~l'^ z 

h^O h V27T 

^(i/--)ie- 2 / 2 dz. 



Since [i is a piecewise linear function, and = liminf 



2-2/ 



it's easy to get 



/OO Z 1 2 
) e~ 2 I 2 dz > //.(y), VyGR. 
-oo n V27T 
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Hence, linij^oo J Q * n' n (X s )ds > Jq fj,'_(X s )ds almost surely. Applying Ito's formula to G n (X t ), 
we have 

[\ n (X s )dW s = G n (X t ) - G n (X ) - l\ n (X s ) 2 ds-\ f\' n (X s )ds. 
Jo Jo 1 Jo 

Passing to the limit as n — > oo yields 

I r l (X s )dW s < G(X t ) - G(X ) - f v(X s ) 2 ds - \ I n'_(Xs)ds, a.s. (2.4) 
Jo Jo 1 Jo 

Finally, combining (|2.4p with (|2.3j) we have 

pit, x, z) 

which yields the desired result (|2.1|) . ■ 

Lemma 2.2 For x, z < c, h > 0, it holds 

P x ( sup (X s - c) < 0\X t = z)< 4^4e G W- G (-)- 3 ^/ 2 (l - e - 2(c "1 (c ' z) ), (2.5) 
0<s<t P{t,x,z) 

and 

P.( sup (X s -c) > 0|X t = z) < l-g p^ + C ~* ) + e ^(^)$( -^ +/t - c+z ), (2.6) 

t<s<t+h v Vft 7 V Vrt 7 

where fi + =sup J/gR /x(y), 3>(y) = J^o ^^ e_z2 ^ 2 ^ z ' anc ^ ^(') are defined in Lemma [2A\ 



Proof, (i) By Anderson [I] or Siegmund [8] (P. 375 ), we have 



"2( c x)(c z) 

Q*( sup (W s - c) < 0) = 1 - e i . 

0<s<t 



Substituting the above into f|2.1|) results in f|2. 5j) . 
(ii) Define a new process 

dX s = // + ds + dW s , X t = z. 
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Then according to the comparison theorem ( Ikeda and Watanabe [5j ), conditioning on Xt = z, 
X s < X s for every s > t almost surely. Therefore 

P x ( sup (X, - c) > 0\X t = z) 

t<s<t+h 

<P( sup (X s -c)>0) 

t<s<t+h 

= P( sup - t) + W, - Wt + z-c) >0) 

t<s<s+h 

= Q ( sup (W s + n + s + z-c)>0) 

0<s<h 



where the last equality follows from Anderson [T]. This yields ()2.6p . ■ 

Theorem 2.3 (Existence of density) Let (Xt)t>o be defined by U.l\) with //(•) satisfying 
tl.3\) and r c defined by J l.^p . T/ien r c /ias a density f c (t,x) satisfying 

n / \ C — X CV„1 3Mt (g-^) 2 ,„ 

/c( *' x) " 7»7i e e "^' (2 ' 7) 



where x < c is the initial state of the process (Xt)t>o- 

Proof. In order to show the existence of the density of r c , by Radon-Nikodym theorem we only 
need to prove that, for every h > 0, 

limP a; (r c G (M + h))/h < oo, 
/k|,0 



where 



P x {T c e{t,t + h))= f P x (T c e{t,t + h)\X t = z)P x (X t edz) 

+ / P*(r c G (t,t + fc)|*t = «)P s (X t G dz) 

-oo Jc-h 1 / 4 



Il+I 2 
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To this end, we estimate Ii and I2 subsequently. First, due to (|2.6j) we have 

Ii= / P x ( sup(X s - c)<0\X t =z)P x ( sup (X s -c)> 0\X t = z)P x (X t £dz) 

J-00 0<s<t t<s<t+h 

< P x ( sup {X s - c) > 0\X t = z)P x (X t G de) 

J-00 t<s<t+h 
c 



Now we show that / e 2fl+( - c ' z) P x (X t £ dz) < 00. Let X s be the solution of the SDE 
J— 00 

dX s = n-ds + dW s , X = x, 

where [i- = infygR fi(y) which is finite by assumption. Then Xt < Xt almost surely for every 
t > 0. Hence, for each A < 0, 

E(e AXt ) < E(e A ^) = E(e Xfl - t+XWt ) < 00. 

By a similar argument, we have for each A > 0, 

E(e XXt ) < E(e x >*+ t+XWt ) < 00. 

Thus, we have 



J— 00 



Moreover, since 



1 - g( ^ + ^ = o( h ), *( -^ h - ^ = o(h), as U0, 



it follows that Ii = o(/i) as /i 1 0, where o(/i) is the little-o notation which means o(h)/h — > as 
0. 
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Next we estimate li- Using Lemma 12. II and Lemma |2.2| we have 

h = [ P«( sup(X,-c)< 0\X t = z)P x ( sup (X s - c)>0\X t = z)P x {X t edz) 

Jc-h 1 / 4 0<s<t t<s<t+h 



< / " q(t, x, z)e G ^- G ^~ ^ (1 - e -2(—)( C -)A) 



dz 



o/ \ i-h x l 4 

= 2(C-Xj / e G( e -y)-G(*)-2^ (y + 0(y 2 )) r i _ $( _^l/ 2 + ^-l/2) 
* Jo L 

+ e 2 ^<&( - n+h 1 ' 2 - yh- 1 / 2 )] q(t, x, c - y)dy, 
where 0(y 2 ) is the big-0 notation. It follows that 



1 



rh 1 / 4 



lim l 2 /h < Klim^- I (y + 0{y 2 )) \l - $( - fi + h 1/2 + y/i~ 1/2 ) 



h\0 



(2.8) 



+ e 2 ^+ y ^ ( - ^+/i 1/2 - y/i~ 1/2 ) 1 g(t, s, c - y)dj/, 

where X = 2 i^LJ!l e G ( c )~ G ( x ) — . To obtain the limit of the right hand side of ()2.8j) . we only 
need to calculate 



1 r hl/4 r 

\\ y[l-*(-n + hV 2 +yh-V 2 ) 



lim — / y 

hio h J 



+ e 2 ^+ y $( - fi+h 1 / 2 - yh- 1 / 2 )] q(t, x,c- y)dy 



= lim / y— 

HO J dh 



l_$(_^ +/l l/2 + ^- 1 / 2 ) 

+ e 2ll + y <S>( - n + h 1/2 - yh- 1 / 2 )] q(t, x,c- y)dy 



rh 1 ' 4 



= lim 





2^2^ L 








+ e 




rh- 1 ' 4 


lim 




hiO 


Jo 



^ + h-'l 2 + yh-^ 2 ) 



( - /i+Zi" 1 / 2 + yh-^ 2 ) 1 g(t, x, c - y)dy 



(M+h 1/2 +;;) 2 . , /9 . 



+ e ^ ( - fi+h 1 ' 2 + y) 



(c-x-yh 1 / 2 ) 2 

e at dy 
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Substituting this into (|2.8p . we have 



l iml2 / h < _^_ e G( C )-G(x)-^ e -^ 
HO 1 ~ V2^t 3 / 2 

Analog to the arguments on Ii, we finally obtain 

1 . , ,nn ^ C — X rlr \ rl , A 3Mt (c-x) 2 

lim-P x (r c E (t,t + h)) < G(c)-G{x)- — e gj- 

Hoh y v^t 3 / 2 

which yields the existence of the density f c (t, x) and its upper bound (|2.T|) . ■ 

After prooving the existence of the density for the first passage time r c of diffusions with 
piecewise linear drift, now we use the method of Darling and Siegert [2] to obtain the distribution 
of r c in terms of its Laplace transform. Recall that p(t, x, y) denotes the transition density of 
the process (X t )t>o and f c (t,x) is the density of r c for the process (X t )t>o starting at Xq = x. 
Their Laplace transforms are defined as 

/•oo 

p(x,y\X)= e~ xt p(t,x,y)dt, 



o 



oo 

Xt 



fc(x\X)= / e- xt f c (t,x)dt. 
Jo 

In [2], the authors study the properties of p and f c under the assumption of the existence 
of f c and give the following results. 

Proposition 2.4 (|2j) (i) p(x,y\\) is in a form of product 

u{x)ui(y) y > x 



P{x,y\X) = < 



(2.£ 

v{x)vi{y) y < x, 



and 



f c (x\X) = ^l, x<c. (2.10) 
u{c) 



(ii) the functions u{x) and v(x) can be chosen as any two linearly independent solutions of the 
differential equation 

with initial condition u(—oo) = v (oo) = 0. 
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Proof. The key points of the proof are based on the following basic observations: 

P(t, x, y) = \ f c (s, x)p(t - s, c, y)ds, for x < c < y; 
Jo 

dp dp 1 d 2 p 

and p{t,x,y) satisfies: — = ^{x) — + 2 a ^~dx 2 ~ ' 

We refer to [21 Theorem 3.1 and 4.1] for the complete argument. ■ 

According to the theory of ordinary differential equation, the explicit solution of equation 
(|2,lip can be given for some special cases. Therefore one feasible method to find the first 
passage time distribution of a general diffusion process is to approximate this process by a 
kind of diffusion processes for which the corresponding equation (|2.11|) admits explicit solutions. 
Based on this idea, we establish the following main theorem which provides an explicit solution 
of (|2.1ip for a diffusion process (X t )t>o defined by (jl.ip and has drift f)1.3|) . 

Theorem 2.5 Let (Xt)t>o be the solution of SDE with fi satisfying $1.3\) . Then for x < c, 

u(x) 



fc(x\X) 



u(c) ' 



(2.12) 



with 



u(x) 



ILj \ Ji ( ~ T~ ' o' ~ ai ( x + ~ ) 2 ) la ^° + ji(bi\x)l ai =o l[ Xi _ uXi ](x), (2.13) 

i=— oo 

and satisfying u(—oo) = 0, where 



1 1 13 

Ji{a, ^,x) = C!^(a, - ; x) + C 2 ^x l l 2 ^{a + -, -; x), 

1 °° (a' 

#(a,-;x) = 1 + ^7^^, (o) fc :=a(a + l)...(a + k-l) 



(a)k x k 

k=l V 2^« 



e~ hx [C 1)ie A ^/ 2 + C 2 ,ie- A ^/ 2 ] , if Aj = 46? + 8A > 

JQh; x) = { e -*Hx [ Clji sin(A,x/2) + C u cos(Aix/2)] , if A 2 = -4fej - 8A > 

e-^{Cx,iX + C 2>i ), if A 2 = Ab 2 + 8A = 0. 

Here the constants Cn, C 2) i are chosen such that u{x) is differentiable. They can be determined 
recursively: if C\i, C 2) i for all i < k are known, then Ci^+i, C2,A:+i are chosen such that 
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]im x ^ Xk u(x) = lim xlXk u(x) and linx^ u{x l_ u x ^ k) = lim^ " ( ^_"|f fc) • Furthermore, they are 
uniquely determined up to a multiplicative constant, which will not change the value of f(x\X). 



Proof. The explicit solutions of the following second-order differential equations are known. 
By 2.1.2.11], the general solution of ODE 

■j-o + a T" + by = 
dx z ax 



is 



exp (— \ax) [Ci exp (|Ax) + C2 exp (— |Az)] 
exp (— \ax) [C\ sin(^Ax) + C2 cos(^Ax)] 
exp {—^ax){C\x + C2) 



if A 2 = a 2 - Ab > 0, 
if A 2 = 4b - a 2 > 0, 
if a 2 = 46, 



where C±, C2 are constants, similarly, by [7, 2.1.2.108], the general solution of ODE 



ld 2 y 



dy 



2d~^ + ( aiX + bl ^dx +b ° y = 0, Ql ^ ° 



is 



J (bo, 7;', -ai(x + — ) 2 ), where 
2 a\' ' 



J(a, b; x) = Ci*(a, b; x) + C 2 x 1 - b ^(a - b + 1, 2 - b; x), 
Ma b . x) = 1 + — 

Combining them together, we get the explicit solution f)2. 13|) to the differential equation 



ld 2 y 



Av 



2dx 2+ ^ X) fx- Xv = {) - 



The result follows by applying Proposition 



Remark 2.6 It is worthwhile to include some asymptotic relations of Vl/(a, b;x) here. 



*(a,6;x) = ^rle x x a - b 



ty(a,b;x) 



r(a) 
r(6) 



x 



if x — )• +00, 



r(6 - a) v 7 
where T(z) = J Q °° t z ~ l e~ t o\t is the gamma function. 



l + Of r-r 



if x 



-00, 
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Remark 2.7 Our approach to prove the above theorem can be used to deal with the first 
passage time of two-sided constant boundaries. Similar to our treatment of Theorem 12.31 we 
can prove that the first passage density for two-sided constant boundaries also exists. Then using 
[2j Theorem 3.2], we can get the explicit expression of the Laplace transform of the corresponding 
first passage time density. Inevitably, the expression in this case will be more complicated. 

3 Approximation of diffusions with general drifts 

In this section we consider the boundary crossing probability for general diffusions. The basic 
idea is to use diffusions with piecewise linear drifts to approximate the diffusions with general 
drifts. The key issue is to estimate the accuracy of this approximation. 

To this end we first establish a rule to control the accuracy of the approximation. Our 
following result, Proposition l3.1l below. may be viewed as an extension of Downes and Borovkov's 
result [H Theorem 4.1] to the situation of diffusion processes with non-differentiable drifts. 
Again, let (X t )t>o be a diffusion process satisfying 

dX t = fi(X t )dt + dW t , X = x. (3.1) 

For each e > 0, let (Xf) t >o be a diffusion process satisfying the SDE: 

dXf = He{X e t )dt + dW t , X £ Q = x. (3.2) 

Assume that 

(HI) fj, e : R — > R is bounded and there exists a positive constant K such that 
\Ve(y) ~ Ms 00 1 < K\y - z\, Vy, z G R, e > 0; 

(H2) su P Ky)-/i e (j/)| <e. 

For c > x, define the first passage time for the process {Xf) t >Q as = inf{i > 0; Xf > c}. 
Further, let /i_ = inf ye n fj,(y) and yu+ = sup J/gR /x(y) denote the lower and upper bound of the 
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function fj,(y) respectively. Obviously assumptions (HI) and (H2) imply that [i- and n + are 
finite. 

Proposition 3.1 Let (Xt)t>o and (Xf)t>o be defined as in \3. 1\) and 113.2]) . Then under (HI) 
and (H2), for every T > 0, it holds 



|P(t c >T)-P(t c £ >T)| 



< 2Te 3KT/2 e G(c) - G{x) 



f T c ~ x {c - x)2 ( \ i 1 \ j 1 , , 

I — r e 2s M-H = ds e + o(e), 

Jo ^s3/2 V IF 1 ^27r(T-s)/ J 



where G(y) = fi(z)dz for some fixed yo. 



Proof. First we have 



&(X t -Xl) = (n{X t )-n £ (Xl))c\t 

= (p(X t ) - LL £ (X t ) + fl E (X t ) - ^ £ (X!))dt. 

By (HI) and (H2), 

\X t -Xf\<et+ [ K\X s -X £ s \ds 
Jo 

and it follows from Gronwall's lemma that 

\X t -X £ \ < ete K \ a.s. Vt > 0. 



Using this inequality, we have 



P( sup (X t - c) < 0) - P( sup {X £ t - c) < 0) 

0<t<T 0<t<T 

< P( sup {X £ t - ete Kt - c) < 0) - P( sup (Xf - c) < 0) 

0<t<T 0<t<T 

< P(0 < sup (XI - c) < eTe^ T ) 

0<t<T 

= / P(r c e g ds)P( sup (XI - c) < eTe^ T |X s £ = c). 

Jo s<t<T 

Since //_ is finite by (HI) and (H2), 

X\ -X e s >^(t-s) + W t -W s , t>s, 



(3.3) 



(3.4) 
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which yields 

P( sup (XI - c) < eTe KT \X £ s = c) 

s<t<T 



(3.5) 



(3.6) 



< P( sup (//_ t + W t ) < sTe KT ) 

0<t<T-s 

= -^_(T-s) + sTe KT , _ 2 „_ eTe KT ^ ( -LL-(T-s)-eTe KT ) 
=:I(T, 8 ,K,n-;e). 

Substituting the above inequality into (j3.4j) and using Theorem 12.31 we obtain 
P( sup (X t - c ) < 0) - P( sup (Xf - c) < 0) 

0<t<T 0<t<T 

where G e {y) = ^ e (z)dz, M £ = inf{/^ + 1^ _(y)}, and 

[i (y) :=liminf fe(z) ~ fe(y) <K 

Further, by the inequality 1 — e _x < x for all x 6 R, we have the following upper bound for the 
term I(T, s, K , (j,- ; e): 

f (r, a, M _ ; e) < ( + 2| M _ |Te T jg. (3.7) 

In order to prove the inverse direction of the inequality, we write 
P( sup (XI - c) < 0) - P( sup (X t - c) < 0) 

0<<<T 0<i<T 

< P( sup (Xf - c) < 0) - P( sup (Xf + ere^ T - c) < 0) 

0<t<T 0<i<T 

< P( sup (XI + eTe^ T - c) < eTe^ T ) - P( sup (Xf + Te^ T - c) < 0) 

0<t<T 0<t<T 
/■T 

= / P(r| E ds)P( sup (Xf - c) < eTe^ T |X s e = c), 

JO s<t<T 

where e is small enough so that c — eTe KT > x as c > x, and c = c — eTe KT . Using the same 
arguments for (I3.5P and (I3.6p . we can obtain the lower bound. Combining this with (I3.7P and 
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noting that M £ > —K/3, we have 

|P( sup (X t - c) < 0) - P( sup (X t - c) < 0)| 

0<t<T 0<t<T 

< 2Te ^Tl2 e G { c)-G {x) [ ^ T _g_ e -^# ^_ | + i/(2vr(r - 8 ))i)ds] e + o{e), 
which concludes the proof. 



Let Xt be the solution of SDE (|3.ip . Assume that there exist positive constants M\, M 2 
such that 

sup \n(y)\ < Mi, sup \u(x) - n(y)\ < M 2 \x - y\. (3.8) 

For each n £ N, we partition the real line (—00, 00) with sub-intervals with endpoints {0, ±-, ±~, . . .}. 
Let // n (y) be the piecewise linear function taking values at xi = ~ (i G Z) to be fi(xi). Thus, 
for each j £ Z, 

max \u(y) - ri n (y)\ < max { max{|^(i/) - fJ,(xi)\, - n(x i+1 )\}} 

y£[xi,Xi +1 \ ye[xi,x i+1 \ ^ 

< M 2 /n. 

(n) 

Further, let X± be the solution of the SDE 

dX f (n) = u n {x[ n) )dt + dW t , X^ n) = x. (3.10) 

Then for each n 6 N, (X^ ) t >o is a diffusion process with piecewise linear drift. For c > 
x, Theorem 12.51 gives the distribution of its first passage time r c (n) := inf{t > 0; xf l) > 
c} in terms of its Laplace transform. As we have stated at the beginning of the paper, the 
boundary crossing probability of the general diffusion process (X t )t>o can be approximated by 
the boundary crossing probability of (X t )i>o- The accuracy of this approximation is given in 
the following theorem. 

Theorem 3.2 Let (Xt)t>o and (X^"^>o be the processes defined above. Then under condition 
( TOP . forc>x,T> 0, it holds 

|P(r c > T) -P(r c ( n ) > T)\ 

< 2TM 2 e 3M * T / 2 e M ^- x ~> \ C " l e~^ f ^ + -=^_) dal - + o(-). ' ' ' 
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Proof. By the definition of fJ, n (y), it is easy to see that 

sup \Hn(y)\ < Mi, \(i n (y) - fi n (z)\ < M 2 \y - z\, y,z£~R. (3.12) 
yen 

It follows that G(c) — G(x) < Mi(c — x). Therefore (|3.1ip follows from Proposition 13.11 and 

(H2D. ■ 
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